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ABSTRACT 

Inhomogeneous reionization acts as a source of arcminute- scale anisotropics in the cosmic mi- 
crowave background (CMB), the most important of which is the kinetic Sunyaev-Zel'dovich 
(kSZ) effect. Observational efforts with the Atacama Cosmology Telescope (ACT) and the 
South Pole Telescope (SPT) are poised to detect this signal for the first time, with pro- 
jected 1 /iK^ -level sensitivity to the dimensionless kSZ power spectrum around a multipole of 
/ = 3000, [A/3ooo]^- Indeed, recent SPT measurements place abound of [A/3000]^ < 2.8 /iK^ 
at 95% C.L., which degrades to [A^aooo]^ < 6 /iK^ if a significant correlation between the 
thermal Sunyaev-Zel'dovich (tSZ) effect and the cosmic infrared background (CIB) is al- 
lowed. To interpret these and upcoming observations, we compute the kSZ signal from a suite 
of ^ 100 reionization models using the publicly available code 21CMFAST. Our physically 
motivated reionization models are parameterized by the ionizing efficiency of high-redshift 
galaxies, the minimum virial temperature of halos capable of hosting stars, and the ionizing 
photon mean free path - a parameterization motivated by previous theoretical studies of reion- 
ization. We predict the contribution of patchy reionization to the / = 3000 kSZ power to be 
l^moo'^? = 1-5-3.5 /iK2. Therefore, even when adopting the lowest estimate in the litera- 
ture for the post-reionization signal of [A[^^qq]^ ^ 2 /iK^, none of our models are consistent 
with the aggressive 2a SPT bound that does not include correlations. This implies that either: 
(i) the early stages of reionization occurred in a much more homogeneous manner than sug- 
gested by the stellar-driven scenarios we explore, such as would be the case if, e.g., very high 
energy X-rays or exotic particles contributed significantly; and/or (ii) that there is a signifi- 
cant correlation between the CIB and the tSZ. The later is perhaps not surprising, as massive 
halos should host both hot gas and star forming galaxies. On the other hand, the conservative 
SPT bound of [A^aooo]^ ^ 6 /iK^ is compatible with all of our models and is on the thresh- 
old of constraining physically motivated reionization models. The largest patchy kSZ signals 
correspond to an extended reionization process, in which the sources of ionizing photons are 
abundant and there are many recombinations (absorptions in sinks). We point out that insights 
into the astrophysics of the early Universe are encoded in both the amplitude and shape of the 
kSZ power spectrum. 

Key words: cosmology: cosmic background radiation - dark ages, reionization, first stars - 
early Universe - diffuse radiation - large scale structure of Universe - early U 



1 INTRODUCTION 

Scattering of cosmic microwave background (CMB) photons off 
of inhomogeneities in the electron density after the recombination 
epoch imprints "secondary" anisotropics in the CMB. These sec- 
ondaries are the dominate contribution to the CMB anisotropics at 
< 5^ scales. The largest and most studied of these secondaries is 
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the thermal Sunyaev-Zel'dovich (tSZ) effect, which results from 
Compton scattering off of ^ 10^ K intracluster gas (Zeldovich & 



Sunyaev 1969 ), and most of the tSZ arises from z < 1 (e.g., |Ko- 
,matsu & Seljak^2002 ). The second most prominent is the kinetic 
Sunyaev-Zel'dovich (kSZ) effect and is the subject of this study. In 
contrast to the tSZ, the kSZ is not sourced by rare regions in the 
Universe but arises because most of the volume has a bulk pecu- 
liar motion with respect to the CMB frame ( Sunyaev & Zeldovich] 
[T98Ql[Ostriker & Vishniac|1986[|Vishniac,1987„Ma & Fry,20021 ). 
These motions impart a Doppler shift to the 10% of CMB pho- 
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tons that scatter after recombination. The kSZ is typically generated 
at higher redshifts than the tSZ, with a large fraction likely being 
sourced by inhomogeneities at 2; > 6, during the epoch of reioniza- 
tion (|Gruzinov & Hu|1998l|Knoxet al. 1998 , Valagea s et al. 12001 



Santos et al.||2003| [Salvaterra et al.||2005> ,McQuinn et al.||2005 



Zahn et al.|2005l ). 

Until recently, the prospects for detecting the kSZ were per- 
ceived by the CMB community as bleak. While the unique fre- 
quency dependence of the tSZ anisotropics should allow the sepa- 
ration of the kSZ from the tSZ, theoretical projections were that the 
tSZ signal would be more than an order of magnitude larger than 
that of kSZ on arcminute scales (e.g., Zahn et al. 2005 ). As such, 
instruments targeting the CMB secondaries would face the diffi- 
cult task of using a limited number of frequency channels 3) 
to isolate the much smaller kSZ signal while also using this fre- 
quency information to isolate the contamination from dusty galax- 
ies and active galactic nuclei ( Huffenberger & Seljak 2005 ). How- 
ever, recent measurements have indicated that the amplitude of the 
tSZ power spectrum may actually be comparable to that of the kSZ 
rather than an order of magnitude larger. This lower value has been 
attributed to the lower as that is preferred by current cosmologi- 
cal measurements relative to the estimates of a few y ears ago (the 
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amplitude of the tSZ power spectrum scales as 

|Seljak|[2002] ), and also to the realization that feedback processes 
strongly suppress the tSZ in clusters (Battaglia et al."2010', 'Shaw| 
|et al.|2010 , Trac et al.|2011| ). As a result, separating the two signals 
is now expected to be an easier task than had been anticipated. 

Both the Atacama Cosmology Telescope (ATC|^ and South 
Pole Telescope (SPT|^ began measuring the CMB at arc-minute 
scales in 2007 and have recently published upper limits on the 
kSZ angular power spectrum (defined in Section 2). At present, 
the greatest challenge for these efforts is separating out several dif- 
ferent terms that contribute to the angular power spectrum: the pri- 
mordial CMB, the tSZ, radio sources, the cosmic infrared back- 
ground (CIB), in addition to cross-power between the tSZ and radio 
sources, and between the tSZ and the CIB. The primordial fluctu- 
ations are important at low multiples and the sources (CIB/dusty 
galaxies and radio sources) are most important at high multipoles. 
The tSZ-CIB cross-correlation is currently the least understood 
component. 

With act's two frequency bands, jPunkley et al.| ( |2011| 
placed a limit of / [27v][Ct^^ + Cf] = 6.8 zb 2.9/11? 
at the angular multipole of / ^ 3000, where Cf^ and C^^^ 
are the tSZ and kSZ angular power spectra, respectively. Using 
SPT's three frequency bands and ignoring the tSZ-dusty galaxy 
correlation, |Reichardt et al.| ( [2011^ ) recently re ported an upper 
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bound /V[27r]Cf^^ < 2.8 /iK^t 95% C.L.. 
(|2011 ) did note that the kSZ constraint is significantly weakened 
to /^/[27r] < 6 /xK^ if this study allows for a significant anti- 

correlation between the tSZ and the cosmic infrared background. 
Analysis of the Planck data ( |Ade et al.|2011| ) suggest that 90% of 
the dust anisotropy power at / ~ 3000 is from z > 2 while simula- 
tions of the tSZ signal predict that the bulk of the signal is sourced 
from lower redshifts ( Battaglia et al.|2011|[T7ac et al.|2 01 1 ). While 
this Planck analysis argues that the tSZ-CIB cross-correlation is 
small, further work is needed to better characterize the redshift dis- 
tribution of these signals as current models do not provide good fits 
to the Planck data ( jShang et al.||2011| . The detection of the kSZ 



signal will likely require a combination of Herschel measurements 
of the dusty galaxies and multi-frequency millimeter observations 
by SPT and ACT. Nevertheless, even the recent conservative bound 
of /^/[27r] C^^^ < 6 fiK^ rules out the most extended reionization 
models in |McQuinn"etarH2005l ). 

The kSZ is comprised of two components: a component 
sourced primarily aft er reionizatio n and referred to as the Ostricker- 
Vishniac (OV) effect jOstriker & Vishniac l986,Ma & Fry 2002p] 
and a component from during reionization referred to as the patchy 
kSZ ( Gruzinov & Hu|1998[|Knox et al. '1998V The OV effect results 
from scattering off of density inhomogeneities with bulk peculiar 
velocities, with the fractional contribution per unit redshift peak- 
ing Sit z ^ 1 and monotonically decreasing with redshift above 
this peak. The OV is sourced by relatively linear scales, and it is 
likely that it can be modeled well enough to isolate the kSZ con- 
tribution from reionization. The patchy kSZ signal results from the 
order unity fluctuations in the ionized fraction that are anticipated 
by the vast majority of reionization models. Reionization also im- 
prints other anisotropics in the CMB, primarily in polarizati on ([Hu| 
[20001 pore et al.|[2507l [Dvorkin & Smith||2509l |Su et al.||2011| l, 
which are not discussed here. However, other than the large-scale 
(I < 10) polarization anisotropics, the kSZ is by far the most de- 
tectable of the reionization-induced anisotropics. 

The patchy kSZ anisotropics were originally calculated for 
several toy models of patchy reionizati on {(Gruzino v & Hu 1 998[ 
[Knox et al.| 19981 [Valageas et al.|2001[ [Santos et al . 2003 ). In the 
last several years, the kSZ from reionization has been calculated 
for more physically motivated models of reionization either using 
numerical simulations or semi-analytic models ( [Zahn et al.|2005[ 
McQuinn e t al.|2005]|Iliev et ar2007). Numerical simulations of 
reionization have not achieved large enough box sizes to both re- 
solve this process and also capture the velocity flows which are 
correlated over hundreds of Mpc and which contribute significantly 
to the amplitude of the kSZ anisotropics (although, see [Iliev et al.| 
[2007 1 for a method to account for flows larger than the box size). 
More fundamentally, the heavy computational requirements of nu- 
merical simulations do not allow for exploration of the large as- 
trophysical parameter space associated with reionization. Because 
of these limitations, recent years have seen the development of 
so called "semi-numerical" simulations of reionization, which in- 
corporate an excursion- set model of reionization yurlanetto et al.| 
[2004] ) with large-scale realizations of the density field to produce 
3D reionization morphology ( Zahn et al.'2005VMesinger & Furlan- 



ettol[2007l ^Geil & Wyithe 2008, Alvarez et al. 2009, Choudhury 



et al.|2 009, Mesinger et al. 2011 , Zahn et al. 2011 ). In this study, 
we use the publicly available 21CMFAST cod^to generate ~ 100 
reionization scenarios in order to explore the parameter space of 
the patchy kSZ signalj^ 

There are many other potential probes of reionization, includ- 
ing: the CMB Thompson scattering optical depth, Te, (e.g., |Ko-| 
[matsu et"ari[2011|); the redshifted 2 1cm line from neutral hydro- 
gen (see, e.g., [Furlanetto et al.|2006 |and references therein); trans- 
mission statistics of quasar (QSO) spectra (e.g., [Wyithe & Loeb| 



^ In some prior studies, "OV" has been used to refer to the linear-theory, 
post-reionization kSZ and the nonlinear theory was referred to as the "non- 
linear kSZ." We do not adopt this terminology here. 
^ http://homepage.sns.it/mesinger/Sim.html 

^ When this work was nearing completion, a complementary study by'Zahn] 



^ http://www.physics.princeton.edu/act/ 
^ http://pole.uchicago.edu/ 



|et al.lpOll) appeared in the literature. This study focused on constraining 
the duration and timing of reionization with recent SPT results. We will 
briefly compare our results to theirs in ^ 
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20041 [Mesinger & Haiman|2004| [20071 [Fan al.|2006[|Bolton & 



Haehnelt 2007a, Maselli et al. 2007 ; Gailera ni et al.|2008(|McGreer 
et al. 2011 , Bolton et al. 2011 ); intergalactic damping wing ab- 
sorption in gamma-ray burst (GRB) spectra (Miralda-Escude 1998; 
Barkana & Lo eb^'20041. 'Tot ani et al.||2006l jMcQuinn et al.„2008l 
Mesinger & Fu rlanetto 2008^; the number density, clustering and 
evolution of Lya emitters (LAEs; Malhotra & Rhoads 2004 , Santos 
et al.l2004UHaiman & Cen|2005 F urlanetto e t al. 2006 , Malhotra & 
Rhoads 2006 , Kashika wa etal.|26o 6 , McQu inn etal.|20 07 , Dijkstra 
et al.|2007|, Mesinger & FurlanettOa2008bi |Iliev et al.|2008,,Dijk-| 



stra et al.|2011| ). However the CMB -based constraints are likely to 
be the only ones in the next decade that are capable of probing the 
highest-redshift part of reionization, z > 10. 

This paper is organized as follows. In ^ we summarize our 
method for estimating the kSZ signal. Our main results are pre- 
sented in ^ namely the kSZ signal from a large range of reion- 
ization scenarios. In ^we compare our results to previous works. 
In ^ we put our results in context with recent and forthcoming 
observations of the kSZ. Finally in ^we present our conclusions. 

Unless stated otherwise, we quote all quantities in comov- 
ing units. We adopt the background cosmological parameters {11a, 
Qm, ^h, n, as. Ho} = {0.73, 0.27, 0.046, 0.96, 0.82, 70 km s"^ 
Mpc~^ }, consistent with the seven-year results of the WMAP satel- 
lite ( [Komatsu et al.|20lT] ). We use the Fourier convention that is the 
standard in modern cosmological studies in which the 27r's appear 
under the dk's. 



as coherent over much larger scales than those of interest, rj/l 
where / ^ 3000 (e.g.. Ma & Fry|[2002| . Here, ry is the confor- 
mal distance from the observer, Vrms (z) is the RMS amplitude of 
the peculiar velocity, and Pxy is the 3D cross power spectrum of 
the overdensity in X with the overdensity in Y, where the relevant 
overdensities are in the ionized fraction, x, and the density, p. Dur- 
ing reionization the patchy term dominates, i.e., the term that arises 
from Pxx ■ In the absence of redshift evolution in the power spec- 
tra and the electron fraction, the integral in equation ^ grows as 
(1 + Zmax)^^^ during matter domination, where Zmax is the max- 
imum integrated redshift. Thus, the same ionization topology at 
z = 10 contributes slightly more to the kSZ anisotropy than this 
topology at z = 7. More importantly, doubling the duration of the 
epoch when ionization structures are of a given scale approximately 
doubles the patchy kSZ signal, cf^^^^^, at that scale. 

As a final remark, we have argued that reionization almost 
certainly occurred in a patchy manner, consisting of ionized and 
neutral regions. Even without any knowledge of Pxx other than 
that reionization was patchy - composed of zeros and ones -, there 
is an integral constraint on the reionization contribution to C^^^ 
because the variance of the ionized fraction (or — 1, where 
Xi is the ionized fraction) must equal f k^dk/{27T^) Pxx- Thus, 
fixing the reionization history, f dl [/^p^^^yj^ ^ single num- 
ber that is independent of the reionization morphology. This con- 
straint suggests that smaller bubbles result in a smaller signal since 
j^patchyj2 _ (characteristic/)-^ - bubble size. However, the 
clustering of bubbles can invalidate this simple interpretation. 



2 ESTIMATING THE KSZ SIGNAL 

The CMB temperature anisotropy in the direction of the line-of- 
sight unit vector, u, can be written as: 



^(vi)^cjT j dzc{dt/dz)e-^^^'^ 



rieU-v, (1) 



where ax is the Thomson scattering cross section, re{z) is the 
Thomson optical depth to redshift z in the direction u, v(u, ^) is 
the peculiar velocity, and ne(u, z) is the electron number density. 

In our investigations, we split up the integral in eq. ([TJ into the 
post-reionization, low-redshift component - the OV signal - and 
the high-redshift, patchy kSZ reionization component. The patchy 
kSZ signal is defined in this study as the total kSZ from z > 5.6, 
even in cases for which reionization completed at redshifts higher 
than z = 5.6. However, the patchy reionization signal typically 
dominates the total kSZ from z > 5.6 and so this somewhat arbi- 
trary definition makes little difference and facilitates the compar- 
ison of models. While the patchy kSZ signal is the focus of this 
paper, we discuss both the patchy and OV components in turn be- 
low. 

We focus on characterizing the kSZ angular power spectrum 
^cmbl<^T(A:)|^, where St is the Fourier transform of 6t 
and Tcmb — 2.73 K is the temperature of the CMB. This statistic 
is the final analysis product of the relevant CMB observations. We 
will use the notation [A^^]^ = /^Cf ^^/(27r). It is intuitive to 
relate in the flat sky approximation Cf^^ to the 3D ionization and 
density power spectra: 



QkSZ 



^ cmb 



/ 



+ 2Pxp{z, 



dv 2 -2 

a c 



/ n2 2 - 
[Z] CTxrie 



(2) 



2.1 kSZ from patchy reionization 

To model the patchy kSZ signal, we use the publicly available 
21CMFAST cod^ This code uses perturbation theory (PT) to gen- 
erate the density and velocity fields, and PT plus the formalism 
of [Furlanetto et al!] p004| ) to generate the ionization fields. We 
summarize these calculations below. For further details on these 
algorithms, see [Mesinger & Furlanetto] f2007| ) and [Mesinger et al.] 
(^11). 

We first create a Gaussian random field in a simulation box 
which is 500 Mpc on a side, sampled onto a 1800^ grid. This field is 
then mapped onto a lower-resolution 450^ Eulerian grid at a given 
redshift using first-order Lagrangian PT ( Zel'D ovichj 1970 1) to com- 
pute what we refer to as the "nonlinear density field" and corre- 
sponding velocity field. The statistical properties of the density and 
velocity fields have been shown to match those from a hydrody- 
namic simulation remarkably well over the range of scales and red- 
shifts relevant for patchy kSZ measurements ( Mesin ger et al.|20lT| |. 

Our ionization fields are generated directly from the non- 
linear density field, with the excursion- set pres cription of [Furlan-[ 
[etto et al.[ ( [2004[ ), sHghtly modified according to [Zahn et al.[ < |2011| l 
(see their 'FFRT' scheme). This algorithm marks regions as ion- 
ized for which the number of ionizing photons produced is greater 
than the number of neutral atoms. Specifically, a simulation cell at 
coordinate x is flagged as ionized if 



/coii(x, 2;,i?) > C ^ 



(3) 



where we have also approximated the cosmological velocity flows 



where C is an ionizing efficiency parameter (described later) and 
/coll is the fraction of mass residing in dark matter halos inside 
a sphere of radius R and mass M = 4/37ri?V. where p = 
p[l + (^NL)i?]. This tabulation includes all halos above a set virial 
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5.60 6.06 6.58 7.15 7.80 8.52 9.33 10.25 11.30 12.51 13.90 15.51 17.41 19.65 22.33 



Figure 1. Slices through our ionization fields with thickness 1.1 Mpc. The top panel corresponds to {C, ^vir, ^mfp} = {32, 10^ K, 60 Mpc}. The middle two 
panels illustrate the impact of temporal evolution in our astrophysical parameters: {C, T^ir, ^mfp} = {32[(l+2;) /ll]^, 10^ K, 60 Mpc}, and {C, T^ir, ^mfp} 
= {32, 10^ K, 60 [7/ (1 + z)]^ Mpc}, from top to bottom. The bottom panel corresponds to an extreme scenario with a small mean free path, {(", Tvir , Rmip} 
= {32, 6.3 X 10^ K, 3 Mpc}. The models have Te= 0.087, 0.086, 0.085, and 0.098, consistent at la with WMAP, and [Af^QQQy]^= 2.4, 3.3, 2.5, and 3.1 /iK^ 
(top to bottom, respectively). The vertical lines that are sometimes evident (such as in the second panel at 2; ~ 8.3) owe to how the snapshots are stacked (see 
text for details). 



temperature threshold, Tvir. We do not explicitly compute indi- 
vidual halo locations (as in Mesi nger & Furlanetto|20"07| ), but in- 
stead use the non-linear density field when computing /coii jZahn 



etal, 



201 Following the excursion-set approach (e.g., [Bond et al. 



1991[ [Lacey & Cole|[l993| ), we iteratively decrease the scale R, 
starting from some maximum value, i^max, which we specify be- 
low as Rmfp- If at any R the criterion in equation ([3| is met, this 
cell is flagged as ionized. We also take into account partially ion- 
ized cells by setting the cell's ionized fraction to C/coii(x, i^ceii) 
at the last filter step for those cells which are not fully ionized. The 
resulting ionization fields are found to agree well with those gener- 
ated with cosmological radiative transfer algorithms on scales rele- 
vant for the kSZ (Zahn et al. 2011 , Mesinger et al.|2011| ). 
We parameterize reionization with 3 free parameters: 

(i) C- the ionizing efficiency ofhigh-redshift galaxies. This quan- 
tity can be defined as C = fescf^N^/(l + rirec), where /esc is 
the fraction of ionizing photons produced by stars that escape into 



the intergalactic medium (IGM), /* is the star formation efficiency, 
iV-^ is the number of ionizing photons per stellar baryon, and rirec 
is the mean number of recombinations per baryon. For reference, 
/^,^ = 0.1, /* = 0.1, = 4000 (appropriate for PopII stars), 
and rirec = 1 yield C = 20. However, the parameters /esc and /* 
are extremely uncertain in high-redshift galaxies (e.g., |Gnedin et al.| 
[20081 [Wise & Cen|2009[|Paardekooper et al.|2011| ). In this work, 
we explore the range C ~ 10-50, which corresponds to what is con- 
sistent with CMB and Lya forest constraints on reionizationQ 



^ Simple models at moderate redshifts find that a linear scaling between 
halo mass and UV luminosity, like the one used here (albeit for lower halo 
masses), provides a good fit to the observed galaxy luminosity function and 
their clustering properties (e.g.. Vale & Ostriker"2006[ |Lee et alT]|2009) . 
Nevertheless, it is possible that the ionizing efficiency of galaxies varies 
significantly with halo mass at high redshifts. Since the dominant ioniz- 
ing population at early times likely corresponds to a narrow range in halo 
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(ii) Tvir : the minimum virial temperature of halos that host the 
reionization sources. Much of the reionization hterature sets this 
to Tvir ~ 10^ K, the threshold temperature for efficient atomic 
coohng - corresponding to a halo mass of Mhaio ^ lO^M© at 
z ^ 10. However, Tvir could be smaller than lO^K: The first 
stars were likely hosted by smaller halos with Mhaio = lO^^^M© 
(e.g., [Haiman et aLl|T996l [Abel et al.|[ 200?, 'Bro mm et ar|[2002l ). 
However, star formation inside such small halos was likely inef- 
ficient (with a handful of stars per halo), and was eventually sup- 
pressed by the heating from reionization itself or other feedback 
proce sses {[Haiman et al. 2000, Ri cotti et al . 2001 , Mes inger et al.| 
|2006[ [Haiman & Bryan 2006 ). Tvir could also have been larger 
than 10 K: Radiative and/or mechanical feedback (e.g., Springd] 
|& Hernquist||2003] ) eventually suppresses star formation inside 
Tvir < 10 K halos, though the details and timing of these pro- 
cesses are not well under stood at hig h redshifts (e .g., Mesinger & 
Dijkstra|2008l|Okamoto et al.|2008||Pawlik & Schaye|2 009ffNev- 
ertheless, it is unlikely that Tvir is much greater than few x 10^ K, 
since these values approximately latch onto the faint end of the ob- 
served galaxy luminosity functions at ^ ~ 6 (e.g., Bouwens et al.| 



p008 , Labbe et al. 2010 , Salvaterra et al.|201 ll|Finlat or et al.|20rT| ) 
Additionally, high Tvir models have difficulty in latching onto the 
slow evolution of the observed emissivity, as inferred from the Lya 
forest at 3 < z < 6 (see Fig.[T2|and associated discussion). Thus, 
we explore the range Tvir = 10^-3 x 10^ K, which encompasses 
all conceivable Tvir. 

(iii) Rmip'- the ionizing photon mean free path within ionized re- 
gions of the IGM. The mean free path to intersect an optically thick 
system (termed Lyman Limit Systems; LLSs) is measured to have 
been Rmip 50 comoving Mpc at z = 6, albeit with order-unity 



uncertainty (e.g., |Storrie-Lombardi et al.||1994 


|Stengler-Larrea 


et al.||1995||Peroux et al.||2003| |Prochaska et a] 


.120091 Songaila 


& Cowief2010V This parameter becomes important in our models 



when the ionized region size R is larger than Rmfp, resulting in 
most of the ionizing photons being absorbed in LLSs rather than 
in the diffuse IGM, thereby retarding reionization. Thus, the small- 
est values for Rmfp have the largest impact in our modelsQin this 
work, we explore the range Rmip = 3-80 Mpc. 

Our parameterization ignores redshift evolution of the above 
astrophysical parameters. Since reionization is expected to have 
been extended, it is in fact likely that there would have been some 
evolution in these parameters. However, we find that only a fac- 
tion of the reionization history contributes the bulk of the patchy 
kSZ signal, making evolution less important. Nevertheless, we also 
investigate a few scenarios that have evolving parameters. 

In order to explore the patchy kSZ signal in this 3 parame- 
ter volume, we have simulated 103 different reionization models. 
The ionization boxes for each model are computed at redshift inter- 
vals of Az = 0.2, spanning the redshift range z = 5.6-22.4. This 



masses (because of the steepness of the halo mass function), a halo mass- 
independent value of C may still be a good approximation. 
^ Note that the combination of mechanical and radiative feedback, which 
by themselves likely have different mass-scalings, could conspire to cre- 
ate a relatively flat (i.e., with weak halo mass dependence) suppression of 
star formation (e.g., Finlat or et al.|20H ). Such an effect most appropriately 
translates to a lower value of C (e.g., by lowering /*). 
^ Note that in our simple model, the mean free path is spatially homo- 
geneous. There are likely inhomogeneities in i^rnfp during reionization, 
though the extent of these inhomogeneities is currently unknown (e.g., 
|Choudhury et al.|2009[|Crociani et al.(20TT||McQui'iiiret al.,201 1^ . 
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Figure 2. OV dimensionless angular power spectrum for instantaneous 
reionization scenarios with a reionization redshift of Zre = 5.6 (solid 
curve), 8 (dashed), and 11 (dotted). We have also plotted the OV power 
spectrum calculated in linear theory for the Zre = 8 case. 

procedure generates 85 ionization boxes for each parameter combi- 
nation, with the entire ?^100-point parameter space being sampled 
in roughly one week on an 8-core desktop computer. We interpo- 
late the density, velocity and ionization fields linearly in cosmic 
time between two consecutive simulation outputs. As is tradition- 
ally done in secondary CMB anisotropy calculations, we rotate the 
boxes when reaching their edges so as to minimize coherent stack- 
ing of the signal. Sample slices through four simulations are shown 
in Fig. [T] We have performed box size and resolution tests with a 
fiducial model and found convergence in the / ^2000-10000 kSZ 
power to ^ 10% (see Appendix A). 



2.2 The post reionization kSZ 

The kSZ from after reionization - which we refer to as the OV 
effect - is a contaminate for studying the patchy kSZ signal. 
Any modeling uncertainty in this post-reionization component will 
bias constraints on the patchy component. At / ^ 3000, the OV 
anisotropy contributed per unit redshift peaks ai z ^ 1 and falls off 
weakly towards higher redshifts (becoming zero when there are no 
longer free electrons around). Figure [2] shows the OV anisotropics 
for three instantaneous reionization scenarios with reionization red- 
shifts of Zre = 5.6, 8, and 11. We find the amplitude of the OV di- 
mensionless (i.e. without units of length) angular power spectrum 
at / = 3000 to be [Ag^oo]^ = 3-4 fiK^. 

These curves were calculated using the pseudo-nonlinear for- 
mula for the OV derived in |Ma & F7y] ( |2002| ) (a more accurate ex- 
pression than in eg. [2]). T his calculation takes the linear theory OV 
calculation of [Vishniacl ( |1987| ) and substitutes a nonlinear theory 
for the density field, here computed from the Peaco ck and Dodds 



density power spectrum ( Peacock & Dodds 1994| ) |^^| At higher 
multipoles than / = 3000, this suppression is more significant. For 
reference, 7y27r// = {4, 7, 14, 20} Mpc at z = {0.5, 1, 3, 10} 



In our calculation, we use Hnear perturbation theory for the velocity 
field, an approximation justified in |Ma & Fry 1 2002 1 In addition, we note 
that the density power spectrum has not been "filtered" in our calculations to 
account for the pressure smoothing of the baryons (as was done in |McQuinii| 
|et al.|2005"] and references therein), but note that this should have Httle effect 
at the multipole of interest, I = 3000, unless there is significant feedback. 
We verified this for the window function measured in Shaw et al. 201 1 e.g. 



their eq. (21), which we find suppresses the power at I 
few percent. 



3000 by only a 
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at / = 3000. The fact that these are pseudo-hnear scales is why the 
Unear theory calculation only differs from the nonlinear calculation 
by a factor of^2 (compare the dashed curves in Fig. [2] which cor- 
respond to the linear and nonlinear OV for Zre = 8). This relatively 
small difference suggests that the modeling uncertainty is less se- 
vere for the OV compared to, for example, the highly nonlinear 
thermal SZ. 

How do our analytic predictions for the OV compare with 
those in numerical simulations? As mentioned in the introduction, 
simulating the kSZ is notoriously difficult (e.g., [Zhang et al.|2004| l, 
and most numerical calculations of the OV should be considered 
lower limits on this signal owing to missing large-scale velocity 
modes. However, some recent simulations have achieved box sizes 
where they should be capturing the large-scale velocity flows. Of 
note, |Trac et al.| ( [20Tl] ) investigated the OV assuming a similar cos- 
mology as is used here and found a signal that ranged between 
[AgXoo]^ = 1-8 - 2.4 /iK^ at / = 3000 for z^e = 10 in the 
four scenarios they investigated. (The lowest value in this range 
was from a simulation with a strong feedback prescription.) A sim- 
ilar range of values was found in |Shaw et al.| {|2011'), who argued 
that baryonic physics/feedback could suppress the kSZ by a maxi- 
mum of 1 /iK^ from our fiducial estimates that do not include such 
effects. We adopt this number as the maximal suppression in the 
OV from our calculations, i.e., [A^^oo]^ > 2 /xK^. 

We henceforth only plot the kSZ power that originates from 
z > 5.6, which is dominated by the patchy component of the kSZ. 
We have argued that the 2; < 5.6 contribution has [A*^^ 



130001 



3000^Q^''3ooo/2/7r ^ 2-3 /xK^ at / = 3000. This range is 
one-third to half of th e conservative < 6 /xK^ bound of SPT 
i Reichardt et al.||2011 i, resulting in the constraint [Afg^'^QQ^]^ = 
3000^"C^^^^/7r < 3 - However, the more aggr essive 

bound of [Aisooo]^ < 3 /iK^ from Reichardt et al.||201l| would 
imply [Af3-tn'<lMK^ 



3 RESULTS 

3.1 Power spectra and physical intuition 

In this section we present the kSZ signal is several of our models. 
Qualitatively a few trends are anticipated. First, early and extended 
reionization scenarios should have a larger kSZ signal. Second, the 
peak in ionization power shifts from small to large scales as reion- 
ization progresses and the HII bubbles grow (e.g., |McQuinn et aT] 
|2007| ). Therefore the kSZ power at higher multipoles probes earlier 
epochs (see Fig.|4|. 

To quantify these trends, we begin by showing the patchy kSZ 
angular power spectra for a sample of reionization models. Fig. [3] 
shows how the kSZ power spectra change as we vary ( (top panel), 
Rmfp (middle panel), and Tvir (bottom panel). The inset in each 
panel shows the associated reionization history. Unless specified 
otherwise, the parameter values for each curve are {C, Tvir, i^mfp} 
= {42, 10^ K, 60Mpc}. 

Increasing the ionizing efficiency parameter, quickens the 
progress of reionization as galaxies are more efficient ionizers of 
the IGM. The differences in the reionization histories (xhi at fixed 
z) increase with time (decreasing redshift). Therefore changing ( 
has the most significant impact on the smallest multipoles (larger 
spatial scales) that are shown, since these multiples are sourced by 
the later stages when ionized structures were large. 

Changing the mean free path of ionizing photons, i?mfp, pri- 
marily impacts the middle and late stages of reionization. Small 




Figure 3. The dimensionless angular power spectrum of the patchy 
kSZ, [A^^^^^^J^, as C i^mfp^ and Tvir are separately varied around 
{C,Tvir, i^mfp} = {42, 10^ K, 60 Mpc}, (top, middle, and bottom pan- 
els, respectively). The inset in each panel shows the associated reionization 
history. The near- vertical gray dashed line corresponds to the primary CMB 
anisotropies (which are understood well enough that their contribution can 
be isolated at / > 2500). 

values of Rmip extend the end of reionization as seen in the in- 
set in the middle panel of Fig.[3](c.f. |Furlanetto & Mesinger|2009{ 
[Alvarez & Abel|2010| ). When the size of the ionized bubbles ap- 
proaches i^mfp, reionization is slowed as ionizing photons become 
increasingly swallowed by photon sinks instead of pushing out the 
edges of the HII regions. Remaining neutral regions must then be 
ionized by local sources. We see that for this fiducial model where 
atomically cooled halos host ionizing sources, the mean free path 
only affects the kSZ signal when Rmfp < 15 Mpc. As Rmfp is de- 
creased (particularly to the smallest values considered), the peak of 
the patchy kSZ power spectra shifts to smaller scales 

Increasing Tvir delays reionization until the halos above this 
virial temperature threshold can form in sufficient abundance (bot- 
tom panel. Fig. [3]). Increasing Tvir also results in a more sudden 
reionization as the number of halos with time grows more quickly 
as the halos become rarer. 

In Fig. |4] we present the kSZ power spectra corresponding to 
the four models shown in Fig. [T] The total signal is represented 
by the solid red curves, with dashed curves showing the contribu- 
tions up to the specified redshift. The top panels in both figures [T] 
and [4] correspond to a non-evolving model with {C, Tvir, i^mfp} = 
{32, 10^ K, 60 Mpc}. The second panels down correspond to the 
same parameters as the top panels, but including redshift evolution 



^ ^ One might worry that i^rnfp = 3 Mpc is only a factor of 3 larger 
than our cell size. To test whether this resolution limit has a quantitative 
impact on our results, we ran a higher resolution, 750"^, run with the same 
astrophysical parameters as the model featured in the bottom panel of Fig. 
[T] and found good convergence (see Appendix A). 
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Figure 4. The kSZ dimensionless angular power spectra for the models 
featured in Fig.^ The inset in each panel shows the evolution of xhi (y- 
axis) with redshift (x-axis). 
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Figure 5. The cumulative fraction of the z > 5.6 kSZ signal at I = 3000 
sourced by redshifts less than z (bottom panel) and corresponding global 
neutral fractions less than xm (top panel). Three model curves are shown 
corresponding to: {C, T^ir, Rmip} = {32, 10^ K, 60 Mpc}, {32, 6 x 10^ 
K, 3 Mpc}, {80, 2 X 10^ K, 60 Mpc}. 
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of the ionizing efficiency, C = 32 [(1 + z)/ll]^. By increasing 
the efficiency in this model at early times (z > 10) and decreas- 
ing it at later times (z < 10), the ionized structures become more 
uniform in size, and reionization is extended. The corresponding 
power spectra at / = 3000 (corresponding to ^ 20Mpc structures at 
high-z) is the largest in this model, [AJ^3^qqq^]^ = 3.3 /iK^, placing 
it on the threshold of current measurements. Such a model could be 
motivated, for example, by an efficient top-heavy IMF in the early 
universe, transitioning to an inefficient, feedback-regulated star for- 
mation regime at later times. We find that models with stronger evo- 
lution (e.g., C = 32 [(1 + 2;)/ll]^) are ruled out by present QSO 
constraints since reionization fails to complete in time (see our dis- 
cussion of constraints derived from QSO spectra below). Models in 
which ^ decreases with redshift result in a more sudden reioniza- 
tion and a smaller kSZ signal. 

The model in the third panels of figures [T] and [4] also assumes 
a constant C and Tvir as in the top panel. However, this model in- 
cludes evolution in Rmip — 60 [7/(1 + z)]^ . This evolution is in 
rough agreement with extrapolations using the measured trend in 
the ionizing-photon mean free path from z ^2-4 (e.g., Proc haska| 
|et al.|2 009 ), though extrapolations into the reionization epoch are 
extremely uncertain. The evolving i?mfp model results in similar 
HII bubble morphology as the non-evolving model in the top panel; 
however, the end stages are extended. 

Finally, the bottom panel corresponds to an extreme scenario, 
where LLSs strongly suppress the growth of HII regions, resulting 
in a very extended reionization. The corresponding parameters are 
{C, Tvir, i?mfp} = {32, 6.3 X 10^ K, 3 Mpc}. Here ionization struc- 
tures are smaller compared with other models, shifting the peak in 
kSZ power to smaller scales. 

We also point out that there is evidence in all four panels that 
the power in the higher /-multipoles is sourced from earlier times 
when the bubbles were smallest (i.e. the blue curves start "peal- 
ing off" of the red ones from high-/ as the redshift upper limit is 
decreased). As mentioned above, this is due to the fact that the ion- 
ization power peaks on increasingly larger scales as reionization 
progresses and the HII bubbles grow. Understandably, this trend 
is least evident in the bottom panel, where HII bubbles are con- 
strained to be small for a longer duration, and the contrast between 
ionization power on large and small scales is smaller. 

As already discussed, our fiducial models do not include red- 
shift evolution of our three astrophysical parameters. However, it is 
interesting that for most of the models shown in Fig.|4] the majority 
of the kSZ signal is imprinted over a relatively narrow redshift or, 
more fundamentally, xhi range. We have already noted in ^that 
doubling the duration of reionization roughly doubles the signal. 
More precisely however, the important quantity is not the duration 
of the entire reionization process, but just of the epoch when ion- 
ized structures were of the angular scale of interest. The timing of 
this epoch is model-dependent, but we see from Fig.|4]that the ion- 
ization morphology had notable power on / ^ 3000 scales during 
the early/middle stages of reionization. 

To further quantify which epochs source the patchy kSZ sig- 
nal, in Fig.|5]we plot the fraction of the total 2; > 5.6 kSZ signal at 
/ = 3000 sourced by redshifts less than z (bottom panel) and global 
neutral fractions less than xhi {top panel). Three model curves are 
shown, corresponding to: {C, Tvir, i^mfp} = {32, 10^ K, 60 Mpc}, 
{32, 6 X 10^ K, 3 Mpc}, {80, 2 x 10^ K, 60 Mpc}. The extremely 
extended reionization model shown with the blue dashed curve has 
kSZ power imprinted smoothly and uniformly throughout reioniza- 
tion, i.e. d[/S^^Q^^Y / dxm ^ 1. In this case, as is qualitatively evi- 
dent from the bottom panel of Fig.[T] the ionization structure retains 
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power on / = 3000 (L ^ 20 Mpc) scales throughout reionization. 
However, more standard models are not sensitive to the later stages 
of reionization, when HII regions grow to be larger than / = 3000. 
The "fiducial" case shown with the solid red curve only has ^20% 
of the / = 3000 reionization signal imprinted in the last half of 
reionization, xm < 0. 50 Models with even larger HII regions, 
such as the one shown with the dashed green curve, are even less 
sensitive to the later stages of reionization. In particular, this model 
has only ^ 16% of the / = 3000 reionization signal imprinted in 
the last half of reionization. The fact that the kSZ signal is sensitive 
mostly to the first half of reionization suggests that our astrophysi- 
cal parameters can be interpreted as the average values during this 
epoch. 

3.2 Signal at / = 3000 

For the remainder of the paper, we focus our attention of the power 
at / ^ 3000, using the statistic [Af^ooo'']^ = 3000^ Cf^lfoo/ (^^)- 
As discussed in the introduction, these multipoles are the most ob- 
servable with upcoming and future small-scale CMB observations. 
At z > 5, / = 3000 roughly corresponds to a physical distance of 

^ r]27T/l = 20 Mpc. This suggests that [Aj^g^QQQ^]^ depends on the 
timing and duration of ionization structure of comparable scale. 

In Fig.[6] we plot [Af^ooo'']^ as a function of {C, T^ir, i?mfp}. 
Values of [Af^^^^^]'^ were interpolated from our 103 21CMFAST 
samples using tesselation-based linear interpolation. Slices through 
the 3D parameter space are shown at C =15, 52; Tvir = 
^q3.2^ 10^-2 K; and i^^fp = 3 Mpc (left panel). The [Af^ooo'']^ = 
3.5, 2.5, 1.5/iK^ contours are plotted in the right panel. 

In our parameter space the signal ranges from [Aj^g^QQQ^]^ ^ 1 
- 4 /jK^. The strongest signal in Fig. [g] corresponds to the regime 
with early, extended reionization scenarios. In other words, high 
values of [Aj^g^QQQ^]^ result when reionization is driven by very 
small galaxies (small Tvir) whose contribution is significant (large 
so as to start reionization early, and with relatively abundant 
LLSs (low i?mfp), so as to extend the duration of the epochs where 
HII regions are comparable in size to / ^ 3000. These scenarios 
would be the easiest to observe, and the first to be ruled out with an 
upper limit on [Al^^^^^]'. 

There are a couple of interesting non-monotonic trends to note 
in Fig. [5] The peak in power occurs at {C, Tvir, Rmfp} ^ {40, 10^ 
K, 10 Mpc}. As the ionizing efficiency from this peak is increased, 
the kSZ power falls somewhat because this change results in reion- 
ization occurring more rapidly. Lower efficiencies delay reioniza- 
tion to later epochs, also resulting in a slightly weaker kSZ signal. 

Additionally, the smaller the mean free path generally the 
stronger the signal, since they produce more extended reionization 
histories. However, the / = 3000 power is a maximum for mod- 
els with i^mfp ~ 10 Mpc because this Rmip is roughly the scale 
probed at this multipole (see left panel Fig. [6]), whereas smaller 
Rmfp values suppress structure on this scale. 

3.2.1 Including constraints from WMAP and QSOs 

Other probes of the reionization epoch can compliment the kSZ sig- 
nal, and restrict the allowed parameter space. Our two most robust 

Note that in this case, ^ 17% of the 2; > 5.6 signal originates post- 
reionization, 5.6 < z < 8.5, and would be more appropriately classified as 
the OV component (see the discussion of our conventions at the beginning 
of §2). 



constraints at present come from the Lyman- series forests detected 
in the spectra of high-redshift QSOs and from the measurement of 
Te by the WMAP satellite. 

The Lyman-series forests in high-redshift QSOs, originally 
detected in the Sloan Digital Sky Survey (SDSS), have been used 
in recent years to claim that reionization must have been over by 
z ^ 6 (e.g., |Fan et al.|2 006 ). Photons emitted by a QSO blueward of 
the Lya line can redshift into the Lyman resonances along the LOS 
and be absorbed by atomic hydrogen (although complete absorp- 
tion does not necessarily indicate a neutral region since this probe 
saturates easily; | Gunn & Peterson I 1965 ).A model-independent, 
maximally conservative upper limit on xhi can be obtained just 
from the fraction of pixels in the Lyman-a and Lyman-/3 forest 
which ar e jointly dark ([Mesinger|2010b| , yielding xhi < 0.25 at 
z = 5.6 (McGre er et al.|ZU l l! McCJreer et al., in preparation). Fig. 
[T] demarcates in gray the surface corresponding to models which 
have xhi = 0.25 at ^ = 5.6. The parameter space on the other 
side of that surface is robustly ruled out with QSO measurements, 
and we leave it blank in this and future plots. As expected, QSO 
constraints rule out models where the sources reside in very rare, 
massive halos such that reionization is late or those with a signifi- 
cant xhi "tail" to lower redshifts owing to a small mean free path. 

The other significant constraint on reionization comes from 
the measurement of the Thompson scattering optical depth to 
CMB photons. The seven-year WMAP (WMAP7) estimate is Te = 
0.088 ± 0.015 ( Komatsu et al.,201 1). Fig. [7] demarcates the 2a 
WMAP7 limits on Te with purple surfaces and the 3a limits with 
yellow surfaces. We consider the portion of our parameter space not 
enclosed by the 3 a contours to be robustly ruled out, and we leave 
it blank in this and future plots. The WMAP7 Te limits mainly serve 
to rule out very early reionization scenarios. These are the models 
where reionization is driven by very small galaxies. Fig.|7]and Fig. 
ji] demonstrate that an upper limit of [A^^^^^]^ ^ 3.5 improves on 
existing 3a WMAP7 constraints. 

Fig. |8] shows additional slices through the reionization param- 
eter space, to focus on the trends noted above. The top row shows 
[AJ^^QQQ^]^ for cuts at Rmfp = 60 and 3 Mpc. Consistent with our 
previous discussion, [Aj^g^QQQ^]^ increases with decreasing Rmfp 
(until very low values) as reionization becomes more prolonged. 

The middle row of panels in Fig.jsjshows [Af^^^Q^]'^ at C = 21 
(left panel) and 41 (right panel). Contours of constant kSZ power 
shift to more massive sources as the ionizing efficiency is increased 
(i.e. the kSZ signal increases with since reionization occurs ear- 
lier). The non-monotonic trend of [A)^3^qqq^]^ with Rmfp is clearly 
evident in these panels. The signal is not sensitive to the value of 
the mean free path when Rmfp Mpc, since these values are 

large enough that they do not impact the reionization morphology 
and the morphology is driven just by the clustering of the sources. 

The bottom row of Fig. [s] shows slices at Tvir = 10^ and 
10^ K. A similar structure of the signal is seen as in the top pan- 
els, with the signal now decreasing with Tvir- High redshift QSOs 
already rule out most of our parameter space at Tvir = 10^ K. 
Although slightly fainter than current direct detection limits (from 
HubbleAVFC3, see e.g., |Salvaterra et al.|201 l||Finlator et al.|20TT] ), 
these sources reside in halos massive enough to be rare in the early 
Universe. If such relatively massive sources drove reionization, it 
would have to occur at late epochs and occur very rapidly not to 
violate existing constraints, resulting in a weak patchy kSZ signal, 
[^fsooo'^r - l^K^ Even so, these scenarios are still ruled out by 
WMAP7 at > 2a. 

When viewed at constant Te, the reionization signal varies by 
G?[AJ^3^QQQ^]^ ^ l/xK^. If the signal can be determined to such ac- 
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Figure 6. [Aj^g^QQ^-^]^ at / = 3000 from inhomogeneous reionization, as a function of our three astrophysical parameters. Left panel: Slices through the 3D 
parameter space are shown at C = 15 and 52; Tvir = 10^-^ and 10^-^ K; and i^mfp = 3 Mpc. Right panel: Cuts through the parameter space are shown with 
[Ar3'ooon'=3.5,2.5,1.5^,K2. 




Figure 7. 3D parameter space demarcating regions that are inconsistent with current QSO and WMAP7 Te constraints (see text). The gray contour corresponds 
to models which have xui = 0.25 at z = 5.6. The purple (yellow) contours demarcate the 2 (3) a WMAP7 limits on Te. 



curacy (see the kSZ would provide complimentary information 
to the mean optical depth statistic. This complementary information 
can be used to determine the duration of reionization, Azre, as we 
shall see below. In our parameter space, the duration of reioniza- 
tion is most sensitive to the mean free path in the range Rmip ^15 
Mpc. Nevertheless, smaller differences of dlAf^^^^^]'^ - 
can be obtained by varying just Tvir and ( at fixed Te (see the top 
row of Fig.|8]). 

3.2. 2 Redshift and duration of reionization 

Here we cast our patchy kSZ signals in terms of two empirical pa- 
rameters, which can be easier to interpret than our astrophysical 
ones: 



(i) Zre'- defined as the redshift at which xhi = 0.5, 

(ii) A^re^ defined as the redshift interval from xhi = 0.75 
XHi = 0.25. 

Note that by definition Azre is a change in the HI fraction by 
Axhi=0.5. The entire process of reionization could be consider- 
ably longer owing to extended early and late epochs (e.g., |Mesinger| 
|2010a| ), but the kSZ signal at / ^ 3000 is less sensitive to these 
stages as demonstrated above. 

Fig. [9] shows contours of fixed Azre =1.5 and 3.5 through 
our astrophysical parameter space. Understandably, rapid reioniza- 
tion scenarios are located in the top, back right corner of this figure: 
bright sources hosted by massive halos with relatively few LLSs to 
impede the progress of reionization. The contrary is true for slow 
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Figure 8. Slices of [Aj^g^QQ^'^]^ through our astrophysical parameter space. Top: Rmip = 60 and 3 Mpc; middle: C = 21 and 41; bottom: 
K. Also shown are contours at fixed Te, corresponding to the prefereed WMAP7 Te value and its 2a limits. 



10^ and 10^ 



reionization models. For clarity of presentation, surfaces of con- 
stant Zre are not shown; however, they are very similar to the con- 
stant Te contours shown in the previous figures. 



Fig.llOlshows contours of [Af^Q^^^]'^ as a function of Zre and 



Azre- Constant Te lines are also depicted. These almost correspond 
to lines of constant Zre, since the reionization histories of our mod- 
els are relatively symmetric around ZreP^ 



A constraint of [Aj^g^QQQ^]^ < 3.2/iK^ would improve on the 
WMAP7 2a limit on Zre in our models. A detection of [A^f^QQ^]^ 
combined with the measurement of Te could place complimen- 
tary constraints on both the epoch and duration of reionization. 
However, as we have already seen, a measurement of [AJ^3^qqq^]^ 
and Te in general does not uniquely determine Azre- For example. 



Note that the curving of the lower 2cr limit on Te as one approaches 
high values of Azre results from the fact that these extreme scenarios (cor- 
responding to late-appearing sources and small mean free paths), have an 



asymmetric extended tail to low xm (see Fig. [3}. In these models the small 
mean free path to ionizing photons creates an effective horizon thus de- 
laying the completion of the final (i.e. overlap) stages of reionization, as 
discussed above ^Furlanetto & Mesinger|2009^ Alvarez & Abel|2Q10 L 
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Figure 9. Contours of Azre =1.5 and 3.5 through our astrophysical pa- 
rameter space {upper and lower contours, respectively). Note that Azre is 
defined as the interval for which 0.25 < xhi < 0.75, which can be signif- 
icantly shorter than the duration over a larger range in xhi- 




Figure 11. Contours of the logarithmic slope of the power spectrum around 
/ ~ 3000. Specifically, the slope is defined as 2\og[/S^^^^^ / /S^27m^ ^ I 
log[3300/2700]. 




Figure 10. Contours of [AP^qqq^] 2 as a function of z^e and Azre- 

Te = 0.088 and [A^fooo^l^ = 2.^^iK^ allow both Azre = 2.3 and 
4.0. This is primarily due to the fact that the peak of the kSZ power 
shifts in angular scale. 

The slope of the patchy kSZ power at / = 3000 can break this 
degeneracy. The slope requires sensitivities that are a few times 
improved over ongoing measurements th at seek to detect the kSZ. 
This logarithmic slope is plotted in Fig. 



11 



From this, we see 
that the longer/earlier reionization scenarios have a positive slope, 
meaning that the peak in kSZ power occurs on scales smaller than 
/ = 3000, which correspond to earlier stages of reionization (see 
the bottom panel of Fig.|4]and associated discussion). A detection 
of the sign of this slope, combined with Te and [AJ^g^QQQ^]^, should 
robustly determine z^e and Azre- Furthermore, we find that positive 
values of the slope come exclusively from models with low values 
of Rmip ^10 Mpc. Therefore detecting the sign of the slope can 
constrain the abundance of absorption systems in the early Uni- 



verse, at redshifts far beyond those of present estimates based on 
QSO spectra. 



4 COMPARISON WITH PRIOR STUDIES 

Several studies have investigated the kSZ signal from inhomoge- 
neous reionization using either radiative transfer simulations of 
reionization or semi-numeric models of this procesj^ 

• |McQuinn et al.| ( [2005 1 investigated the kSZ from reionization 
using an analytic framework also based on the Furlanetto et aL| 
(|2004) model, by varying the parameter furlanetto et al. ( 2006| ) 
generalized their calculation to include models where the ionizing 
efficiency depends on halo mass. This analytic method is more ap- 
proximate than the ca lculations presented h ere. However, values of 



patchy -| 2 
Z3000 



found in McQuinn et al. 



( 2005 ) are similar to what we 



[A 

find in our analogous models (although, their functional form for 
[AJ^3^QQQ^]^ is slightly more peaked). However, the mass-dependent 
C-models in |Furlanetto et al.| ( |20 06) exhibit a stronger /-dependence 
than we find in our Tvir models, which should mimic the same ef- 
fect. These differences likely owe to the approximations made in 
their analytic derivations. 

• |Zahn et aL| ( |2005| ) calculated the kSZ from two reionization 
models with a fixed and temporally evolving C, parametrization us- 
ing the |Zahn et al.| ( [2Q07] ) semi-numeric implementation of reion- 
ization on top of a 100 Mpc/h SPH simulation. This box size is 
insufficient to capture the largest-scale velocity modes and, thus, 
should underestimate the kSZ signal. Furthermore, it is difficult to 
compare our calculations with theirs owing to their different reion- 
ization historiesp^Nevertheless, we note that their results qualita- 
tively agreed with those of |McQuinn et al.|p005| ) and, by exten- 
sion, those here. 



-"^^ This figure was created with a wider spectral smoothing kernel (A/ = 
600, instead of = 300 used throughout the rest of the paper). This was 
done to smooth over small spectral features present in our realization. We 
find that, unlike [Aj^g^QQ^-^]^, the spectral slope as defined above is sensitive 
to such cosmic variance effects. 



Most previous studies used a cosmology with erg = 0.9, which will 
result in a slightly higher amplitude since the patchy kSZ is proportional to 
v'^ oc cr|. (78 primarily affects the timing of reionization, and not the HIT 
morphology (McQuinn et al. 2007 1. 

The Zahn et al. ( 2007 1 constant C model has a smaller A^^re than ours 
owing to their implementation. 
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• |Iliev et aL] ( |2007| ) used full radiative transfer simulations of 
reionization to calculate the kSZ from a 100 Mpc/h simulation sup- 
plemented with a scheme to account for missing large-scale veloc- 
ity flows. Their power spectra peak at / ^ 3000, similar to our 
constant C models, but have 2 — 3 times higher amplitudes. The 
likely cause of this discrepancy is the presence of a more uniform 
bubble size in the |Iliev et al.| ( [2007] ) models. This is in contrast 
with our models and also with other radiative transfer simulations 
of reionization (e.g., |McQuinn et aL|2007l|Trac & Cen|2007| ). 

• When our work was nearing completion, a similar study was 
made public, focusing on reionization timing and duration con- 
straints implied by the two year SPT data (Zahn et al. 2011). They 
used a reionization model similar to ours, and although detailed 
comparisons are difficult, their power spectra roughly agree with 
ours in both shape and amplitude. However, they do not explore the 
same (C, Tvir, i?mfp) parameter space we do. Instead, they only fo- 
cus on the efficiency parameter, encoding the sharpness of reion- 
ization through an "acceleration" parameter. This effectively allows 
the ionized fraction to rise much more steeply, Xi oc (Co/coii)^, 
where they vary Co and /3 in order to sample a wide range of reion- 
ization histories. Because of this exponential dependence on the 
collapse fraction, they are able to generate much sharper reioniza- 
tion histories, useful for their Markov chain Monte Carlo (MCMC) 
analysis. For example, halving the duration of reionization (and 
roughly the kSZ power) in a fiducial lO^K, Co - 20 
model requires /3 ^ 2.5. Although such a parametrization allows 
for a larger variety of kSZ signals, which can fit the aggressive SPT 
bounds (see below), we caution that it is not physically motivated. 
Any reasonable feedback mechanism delays reionization instead of 
accelerating it (e.g., Thoul & Weinberg 1996 , Springel & HernqmsF 
[2003 , Pawlik & Schaye 200 9).p^Our work here is complementary 
to theirs, as we predict the kSZ signal for a wide range of realis- 
tic reionization models. We also provide physical insights into the 
signal, for example showing that the shape of the patchy kSZ sig- 
nal can be significantly affected by Rmip\ thus [AJ^3^qqq^]^ and Te 
do not uniquely determine the redshift and duration of reioniza- 
tion. We also demonstrate that the bulk of the patchy signal for 
most models is sourced by the early-to-middle stages of reioniza- 
tion. Therefore any kSZ-based inferences about the later stages of 
reionization are indirect and very model dependent. 



5 CONTEXT WITH RECENT AND FORTHCOMING 
OBSERVATIONS 



Reichardt et al. 



(2011 



placed a limit of [A^; 



< 



Recently, 

2.8 /iK^ at 95% C.L. using two year SPT measurements. This ag- 
gressive bound is valid under the assumption that there is no CIB- 
tSZ correlation. Allowing for such a correlation, the constraint de- 
grades to [Aisooo]^ < 6 iiK^ at 95% C.L.. This is also in agree- 
ment with the l yr data from the ACT , which implies [A^sooo]^ = 
6.8 =b 2.9 /xK^ jPunkley et al.|201l| ). 

Despite the wide parameter space of models we explore, we 
predict a fairly narrow range for the patchy reionization component 



The very first, molecularly cooled galaxies could experience posi- 
tive/accelerating feedback due to the enhanced free-electron fraction from 
ionizing photons (e.g., |Oh & Haiman||2002) or hydrodynamical shocks 
( [Shapiro & Kang' 19871, which catalyzes the formation of molecular hy- 
drogen. Such feedback would be mild, confined to fossil HII regions and 
edges of ionizat ion fronts (e.g.,|R icotti et al. 2002 Kuhl en & Madau|2005) , 
and short-lived jHaiman et al.|2000^^Mesinger et al.^2006^ 



of the kSZ power, [APfooo^]^ ^ 1.5-3.5 /iK^ When combined 
wit h est imates fo r the hom ogeneous ionization (OV) component 

this 



in ^2.2 



(see also IShawet al. 2011), [A^^^ 



2-3 iiK' 

work implies that the total kSZ power is [A^sooo]^ ~ 3.5-6.5 /xK^. 
Therefore the recent conservative SPT bound is on the cusp of con- 
straining viable reionization scenarios. 

However, the tighter SPT bound, [A^sooo]^ < 2.8 /xK^ at 
95% C.L., is incompatible with all of our reionization models. 
Smaller patchy reionization signals than our predicted range, i.e. 



chy]2 



< l/xK - as would accommodate this SPT constraint 



- require very rapid and late reionization histories (i.e., spanning 
z ~ 6 — 7). In these, ionizing sources would reside in extremely 
rare, massive halos (Tvir > 10^ K) such that their abundance grows 
rapidly with time. However, such models are difficult to accommo- 
date with current constraints: (i) transmission in the Lyman forests 
of high-z QSOs provides a low-z limit to reionization so that reion- 
ization cannot occur too late (e.g., [McGreer et aL]|2011| ); (ii) the 
redshift evolution of the transmission in the Lyman forests favors a 
slow evolution in the ionizing emissivity and, thus, a more extended 
reionization process ('Mir alda-Escude| [2003 1 [Bolton & Haehnelt| 
2007b ); (iii) the luminosity function of high-redshift galaxies (e.g., 
Bouw ens et al.||2 008 ) supports a sizable contribution to reioniza- 
tion by galaxies fainter and more abundant than those observecf^ 
To illustrate points (i) and (ii) in greater detail, we plot the red- 
shift evolution of the ionizing photon emissivity and xhi in Fig. 
[T2I Increasing Tvir steepens the redshift evolution of the emissiv- 
ity. Therefore, even though one can increase the ionizing efficiency 
to compensate for the paucity of sources (appealing to larger than 
expected values of /* /esc or more exotic, top heavy stellar popula- 
tions), the evolution of the emissivity becomes too steep to match 
the data at 3 < 2; < 6, in the absence of a rapid counter-evolution 
of C or rapid feedback mechanisms. This is illustrated by the dashed 
red curve in Fig. 
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corresponding to Tvir = lO^K, C = 300. 
If the more aggressive SPT bound is confirmed, more exotic 
reionization scenarios that reduce the level of patchiness may be- 
come favorable. The most plausible candidates would be very high- 
energy X-rays with long mean free paths, which would significantly 
"pre-ionize" the IGM (we show above that the early to mid stages 
of reionization are the most relevant for the / = 3000 signal). A 
homogeneous ionization floor of ampli tude Xj suppresses the kSZ 

(see [Visbal & Loeb|2011 



power by the factor ' 



(1 



Xi) 



for a 



simple model of X-ray reionization and the kSZ). For example, if 
X-rays homogeneously pre-ionize the IGM to = 0.3, the patchy 
signal from our models would decrease to [Aj^g^QQQ^]^ < 0.7- 
1.7 /iK^ with some models accommodating the aggressive SPT 
bound. The patchy epoch of reionization would also occur quicker 
given a homogeneous head-start, further decreasing the kSZ signal. 
However, a homogeneous pre-ionization is also an extreme sim- 
plification, as X-rays would still imprint some inhomogeneities on 
angular scales of / = 3000 (^ 20 Mpc at high redshifts), given that 
the majority of X-ray photons have mean free paths much less than 
the Hubble length even for hard QSO spectra. 

In any case, since late-forming massive halos host both hot 
gas and star forming galaxies, it is likely that the CIB and the tSZ 
are correlated (e.g., |Sehgal et al.|2010| ). This would inv alidate the 
aggressive [A^sooo]^ < 2.8 /xK^ constraint proposed by Reichardt 



[et al.[ (|2011j), implying that the more conservative SPT bound of 



Note that the observed galaxies roughly co rrespond to the largest values 
of Tvir ~ few X 10^ K used in this study (e.g., [Labbe et al.|2010|[Salvaterra[ 
let al.|2011||Finlator et al.|2011). 



© 0000 RAS, MNRAS 000, 000-000 



kSZ Signal from Patchy Reionization 1 3 



c 
o 

>^ 
U 

CO 

a 

d 
o 



10 




1 

0.5 




— I — I — I — I — I — I — I — I — I — I — I — I — |— 
A Bolton&:Haehnelt(2007> 
□ McQuinn + (2011) 
T^.^=2xl04K, <-=20 
T^.^=3xl05K, <-=50 
T,.,= 106K, t=300 




•McGreer+(2011): 

I I I ■ I I I ■ I 



10 



12 



14 



Figure 12. The redshift evolution of the ionizing photon emissivity {top 
panel) and xhi (bottom panel). Emissivity constraints inferred from the 
Lya forest are shown as squares (McQu mn et aLpGll} and triangles 
( [Bolton & Haehnelt|2 007b ) (offset by Az = 0.1 to improve legibility); 
the maximally conservative neutral fraction constraints from the dark frac- 
tion in the Lyo; and Ly/3 forests (McGreer et al. 2011) are shown in the 
bottom panel. We also show the analytic predictions, A^ion = C(l + 
^rec) df coll /dt, for three models, assuming rirec = and i^mfp = cxo. 
The solid blue curve corresponds to "standard" parameter choices of Tvir = 
2 X lO'^K, ( = 20, and is in agreement with both constraints. As Tvir is 
increased, reionization occurs more rapidly and at later times. To fit the up- 
per limits on xm, this increase in Tvir must be accompanied by an increase 
in This is illustrated by the green dot-dashed curve, which corresponds 
to the corner of our parameter space, Tvir = 3 x lO^K, = 50. In- 
creasing Tvir also steepens the redshift evolution of the emissivity. There- 
fore, even though one can increase the ionizing efficiency to compensate 
for the paucity of sources, the evolution of the emissivity becomes too steep 
to match the data at 3 < z < 6, in the absence of powerful and rapid 
feedback mechanisms or an equally rapid evolution in C post reionization. 
This is illustrated by the dashed red curve, corresponding to Tvir = lO^K, 
C = 300. 

[Aisooo]^ < 6 /jK^ is more plausible. This limit is consistent with 
all of our reionization scenarios, and is on the cusp of the prob- 
ing early/extended reionization models. A 1 /jK^ constraint on 
the signal is projected for the coming year ( [Reichardt et al.|20Tl] ), 
which should result in the first detection of the kSZ. 



6 CONCLUSIONS 

Observational efforts such as ACT and SPT are beginning to set in- 
teresting constraints on the kSZ power spectrum at / ^ 3000. This 
is highli ghted by the re cent SPT constraint of [A^sooo]^ < 2.8 /iK^ 
( |Reichardt et al.|20TT < 6 /xK^ when allowing for correlations), 
as well as the ACT m easurement of [A^sooo]^ = 6.8 it 2.9 /xK^ 
( [Punkley et al.|20lT| ). A 1 /iK^ constraint on the signal is pro- 
jected for the coming year ( [Reichardt et al.|[20Tl ), which should 
result in the first detection of the kSZ. Second generation instru- 
ments, such as ACTPol and STPPol, will further improve the kSZ 
measurements in the coming years. 

This study calculated the patchy kSZ signal for a suite of 
100 inhomogeneous reionization scenarios, which were gener- 



ated using the publicly-available code 21CMFAST. In these mod- 
els, the contribution of inhomogeneous reionization to the CMB 
angular power spectrum at / = 3000 spans only a factor of four, 
[AJ^3^QQQ^]^ ^ 1-4 /iK^, despite the large volume of parameter 
space covered by the reionization models. The allowed range nar- 
rows further to [Aj^g^QQQ^]^ ^ 1.5-3.5 /jK^ when including cur- 
rent constraints on the reionization history from WMAP and high- 
redshift quasars. We find that for standard models, the bulk of the 
signal is imprinted in the early to middle stages of reionization, 
when the (clustered) ionized structures have angular scales compa- 
rable to / ^ 3000. 

The range [Af^^QQ^]'^ ^ 1.5-3.5 /iK^ combined with esti- 
mates for the post reionization kSZ of [A^^qo]^ ~ ^-3 /xK^, 
results in a total kSZ signal of [A^sooo]^ =3.5-6.5 /jK^. Thus, 
the recent conservative SPT bound of [A^sooo]^ < 6 /xK^, which 
assumes a free, scale-independent CIB-tSZ correlation, is on the 
boarder of constraining reionization, even assuming conservatively 
low values for [A^^oo]^- The strongest and easiest to detect patchy 
kSZ signals correspond to early and extended reionization scenar- 
ios. In such models, the sources of ionizing photons are abundant 
at early times and/or there are many recombinations (absorptions in 
sinks) so that reionization was maximally extended. Furthermore, 
with the assumption of a slope (but not amplitude) in the CIB-tSZ 
correlation, Zahn et alj ( |2011| ) placed an intermediate upper limit 
of [Aisooo]^ < If such a slope is credible, this constraint is 

sufficient to rule out our reionization models driven by molecularly- 
cooled halos (with Tvir ^10^ K), and place limits of Zre ^10, 
Azre < 2 (see Fig.[Tof. ^ ^ 

On the other hand, the tightest SPT bound < 2.8 /xK^ at 95% 
C.L., which assumes no CIB-tSZ correlation, is incompatible with 
all of our reionization scenarios. Although later and more rapid 
reionization models might be in agreement with this constraint, it 
is difficult to push our models to such extremes since these would 
conflict with other observations of the high-redshift Lyman alpha 
forest and galaxy luminosity function (LP). This implies that either: 

(i) the early to middle stages of reionization occurred in a much 
more homogeneous manner than suggested by the stellar-driven 
scenarios we explore, such as would be the case if, e.g., very high 
energy X-rays or exotic particles contributed significantly; and/or 

(ii) that there is a significant correlation between the CIB and the 
tSZ, which invalidates these bounds. Upcoming combined analyses 
of microwave and Herschel far infrared data will likely test the later 
possibility. 

We also find that shape of the patchy kSZ power spectrum 
encodes astrophysical information. An early and extended reion- 
ization in which bubble growth is retarded by recombinations re- 
sults in a kSZ power spectrum which peaks on smaller scales. The 
sign of the slope of the patchy kSZ power spectrum at / ^ 3000 
is useful for discriminating between extended/early and rapid/late 
reionization processes, which have the same [AJ^g^QQQ^]^. In other 
words, the values of Te and [AJ^3^qqq^]^ do not uniquely determine 
the redshift and duration of reionization (c.f., [Zahn et aL]|2011| l. 
However, the sign of the kSZ slope at / ~ 3000, combined with 
Te and [Aj^g^QQQ^]^, is sufficient to determine the redshift and dura- 
tion of reionization. Furthermore, we find that positive values of the 
slope come exclusively from models with low values of Rmip ^10 
Mpc. Therefore detecting the sign of the slope can constrain the 
abundance of absorption systems in the early Universe, at redshifts 
far beyond those of present estimates based on QSO spectra. 

We anticipate that these measurements will be complemented 
by analyses of the Planck EE signal, which is also sensitive to the 
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duration of reionization, and potentially by analysis of the contri- 
bution of patchy reionization to the CMB polarization power spec- 
trum. Reionization probes not based on the CMB will also yield 
complementary information, using 21cm tomography, QSO spec- 
tra and wide-field LAE studies. These tools will soon reveal the 
details of the dawn of cosmic structure. 
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Figure Al. The patchy kSZ power spectra obtained from different simu- 
lations. The (thick) solid blue curve corresponds to the same box size and 
resolution as used in this work: L = 500 Mpc on a side, sampled onto a grid 
with dimensions given by dim=450. The dotted green curve corresponds to 
a larger box, with comparable resolution: L = 750 Mpc, dim=750. The 
dashed magenta curve has the same box size as our fiducial choice, but with 
higher resolution, dim=800.The thin solid blue curves show power spectra 
computed along two basis vectors without rotating the simulation box, for 
our fiducial choices. Finally, the dot-dashed yellow curve corresponds to the 
signal computed without PT, just assuming linear evolution of the cosmo- 
logical fields, such as has been done in previous semi-analytic studies. All 
curves assume {C, T^iv, Rmip} = {32, 10^ K, 30 Mpc}. 
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Figure A2. The kSZ signal corresponding to {C, T^ir, ^mfp} = {32, 6.3 x 
10^ K, 3 Mpc} from our fiducial simulation (solid blue curve and a higher 
resolution, dim=750 simulation (dotted red curve). 

resolution, while the dashed magenta curve has the same box size 
as our fiducial choice, but with higher resolution. Comparing the 
blue curve to the magenta and green ones, suggests that our fidu- 
cial choices result in a slight overestimate of kSZ power spectrum 
for this random seed, but have converged to within 10% over the 
relevant multipoles. 

The thin solid blue curves in Fig. |Al| show power spectra com- 
puted along two basis vectors without rotating the simulation box. 
From their disparity, one can see that the variance of the velocity 
fields on scales larger than our L = 500 Mpc scales is 10% (Fig 7, 
[Iliev et al.|2007| ), requiring box rotations to avoid coherently stack- 
ing the signal. Finally, the dot-dashed yellow curve corresponds to 
the signal computed without PT, just assuming linear evolution of 
the cosmological fields, such as has been done in previous semi- 
analytic studies (McQuinn et al.|2005] l. Ignoring non-linearity re- 
sults in a 30% underestimate of the I ~ 2000-10,000 kSZ power 
in this model. 

One might also worry that our fiducial resolution, with a cell 
size of 1.1 Mpc, is insufficient to accurately capture the kSZ signal 
from our low-i?mfp models. To test whether this resolution limit 
has a quantitative impact on our results, we ran a higher resolution, 
dim=750, run with the same astrophysical parameters as the model 
shown in the bottom panel of Fig.[T] {C, Tvir, i?mfp} = {32, 6. 3 x 
10^ K, 3 Mpc}. We compare their kSZ power spectra in Fig. 
and find that the signal has converged. 
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APPENDIX A: CONVERGENCE TESTS 

Here we briefly quantify the convergence of our numerical results. 
In Fig. |A1| we plot the patchy kSZ power spectra correspond- 
ing to the same astrophysical parameter choices, {C, Tvir, i?mfp} 
= {32, 10^ K, 30 Mpc}, but with different simulation box sizes, 
L, and grid dimensions, dim. The (thick) solid blue curve corre- 
sponds to the same box size and resolution used in this work. The 
dotted green curve corresponds to a larger box, with comparable 
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